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I. INTRODUCTION 



As weapons systems grow in complexity, the cost in time and 
money of developing and testing advanced weapons is so great that mili^ 
tary planners must insure that testing procedures for performance and 
tactics are formulated so as to obtain reliable data for each test. 

Thus the development of valid and reliable models for test design is a 
matter of great importance. Model manipulation is a far cheaper method 
for developing good system design than actual t rial-and-error tests. 

This paper presents a stochastic approach to a combat model of tanks vs. 
tanks with emphasize on the following: 

1. A scheme for development of trial validation criteria. 

2. A useful procedure that can be expanded to a Semi- 
Markov model for combat (land, sea or air combat). 

3. Approachs to data presentation. 

4. Insight into data analysis and data reduction. 

We let a point P, with coordinates (k,£) represent the number 
of casulties of two opposing forces. Thus the coordinates are integers 
0 < k < N and 0 < £ < M, where N is the size (number of units) of 
the friendly force and M is the size of the threat force. It is con- 
venient to normalize this rectangle to the unit square, letting x = k/n, 
y = jt/ra (see figure 1). The trial events are thus points on the unit 
square. Given a point (x,y) for any given weapon system there is some 
probability of transition to some other point (x^,y^), as the result of 

combat. The trial (battle) is over when either force becomes too small, 
that is, our point reaches some predetermined barrier on the top or 
on the right side of the square. It is our belief that if, in any real 
experiment (mock battle), one force quickly overwhelms the other, the 
experiment is a failure, in that it did not really test the system. 
Graphically the point (x,y) would move nearly vertically or nearly 
horizontally across the unit square. Thus we wish to study how to 
design the ' experiment so that the point (x,y) will, with some reason- 
able certainty, remain in a certain size cone about the 45° diagonal. 

In describing our procedure, semi-Markov process and random 
walk terminology are used throughout. We are concerned with the develop- 
ment of the model and its ability to present data. In particular we 
will make certain simplifying assumptions concerning the transition V 

probability matrix for illustrative purposes. 

The technical portion of this paper seems new and provides a 
.discrete state two-dimensional stochastic model- of land, sea or air- 
combat. 



- 1 - 



II THE MODEL 



We assume that the duration of a trial is a random variable 
T and AT denotes the time interval used to update changes in systems* 
state (friendly and/or enemy kills). The realization of times where 
changes in system states are being updated are expressed in terms of: 

0, AT, 2 (AT) , 3 (AT), 4(AT), ... 

AT will be selected based on the following factors: 

1. AT cannot be selected too "small** due to limita- 
tions of the instrument used. 

2. AT cannot be selected so "large** as to permit more 
than one system to be attrited from each force. 



After AT is selected, the particle (representing the systems* 
state) may change its position or stay put as time progresses depending 
upon the results of some corresponding exchange of fire, but once it 
reaches a specified point on the lines x = b and/or y = b (where 
50% < b < 70%), the system or systems are considered to be saturated 
and the problem (engagement-trial) therefore terminated. This termina- 
tion process compares with actual battle conditions where a specified 
force is considered noneffective when a percentage (50-70%) of casualties 
are sustained. (See Figure 1) 



y-axis 




( 1 , 1 ) 



FIGURE 1. 



x-axis 



At the commencement of the engagement (trial), time t, the 
system (friendly and aggressor) is in state (x,y). At the conclusion 
of a AT interval the position of the particle, say (x,y) moves to a 
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position denoted by S^, S^, and as outcomes of the AT time 



interval when the particle moves from (x,y) to , 1_ 

y M 

(* + £- y + S 



X + N> y ’ 



respectively (see Figure 2). can be used to further 



denote the status of the particle-'given no change in the state or incre- 
ment of the particle. Based on the hypotheses used in this model, all the 



(o.oy- 
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FIGURE 2. 



possible positions of the particle are in the unit square and 
outcomes are possible other than the above four, i.e., S q , S^, 

However, instead of using the unit square as the sample space, 
easier to return to the rectangle ft: 



fl = { (x, y) ; 

Further, we define the sets: 

A = { (x,y) : 

B = { (x,y) : 

C = { (N ,M) } 
I' = ( (x,y) : 



x = 0,1,2,. ..,N; y = 0,1,2,. 

0 < x < N, y = M} 
x=N, 0 < y < M} 

0<x<N, 0<y<M} 



no 




it is 



. . ,M} 
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We shall call the totality of points in A or in B or in C the 
absorbing barriers, while all the points in I are called interior 
points. s 

For each point (x,y) in I, we assume we are given the 
probabilities: 

P r* S y* = P» P > 0 

= q, q > 0 

P {S }= r and P{S} = s, r, s^O 
r xy r o 

such that p+q+r+s=l and they are independent of T. For sim- 
plicity assume also they are independent of (x,y). In the actual experi- 
ment, these probabilities are functions of location (x,y). That is to 
say, the transition probabilities depend upon the remaining forces left 
in the battle zone. If one is interested in developing a stochastic 
model for a combat experiment, the probabilities p, q, r, and s will 
generally be functions of time, locations of the particle, and remaining 
units left in the battle zone. 



For each point in the absorbing barriers, we assume that: 



P {S } = 1, 

r o 

P {S } = P {S } = P {S } = 0. 

r x r y r xy 



Denote P(x,y|t) the probability that the particle for the first time 
reaches the point (x,y) at time t and denote P(x,y|t) the probability 
that the particle is in state (x,y) at time t, and let x(x,y) be 
the first passage-time through the point (x,y) namely, the number of 
(AT)'s to be performed until the particle for the first time reaches the 
point (x,y). Then, for any assigned values of N and M, x(x,y) may 
be considered as a random variable, the realization of which equals 
one of the following values, i.e., 



where: 



•t , t + 1, t +2, 
o’ o o 



t + v, . . . (0 ^ v < “) 
o 



Denote: 



t = max (x,y) . 

K (x,y,t) = {k; max(x + y - t, 0) S k < min(x,y)} 
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Denoting by P (i,j,k,£) the probability function of the multinomial 

t 

distribution (p + q + r + s) , the following well-known identity holds: 
P t (i,j,k,i) = TTjr^j p 1 q j r k s l ; a+j+k+i = t) 
and to any interior point (x, y) , i.e., (x,y)el , 



P(x,y|t) = l 



t J x-k y-k k t+k-x-y 

keK (s-k) ! (y-k) !k! (t+k-x-y) ! ^ P 



for t=t,t +1, t +2, where t = max(x,y). 

o o o o J 

For (x,y) and t such that t < max(x,y), it is clear that: 
p (x,y|t) = 0 

For the first-passage probability p(x>y|t) we have the follow- 
ing equations: 

P (0,0 |0) = 1 , 

p (x,y 1 1) =P(x,y|t) - sP (x,y| t-1) . 

If (x,y)el, we should have: 

P(x,y|t) - I (x+lzhii -y-k V- 



keK(x,y 



v M (x-k)! (y-k)! k‘ q P r 
> y 9 L / 



t-i 

x+y-k-l 



t+k-x-y 



provided t S t . 



Proofs of all formulas and optimal mixes of weapon systems used 
for both the friendly force and the threat force will be treated in a 
separate paper. 

Denote by <(> t (u,v) the generating functions of: 

A 

P (x,y 1 1) , namely 
t t 

(u,v) = l l P (x,y 1 1) uV 
x=0 y=0 



Denote $> as the generating function of cf> ^ , then: 
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( 1 ) 



( 2 ) 



co 

4>(u,v,0) = Z 4> t O 11 

t=o 

“ I ® t I I P(x,y|t) U x v^ 
t=o x=o y=o ' 



,-1 



= [1- (qu + pv + ruv + s) ] 

The generating functions of P(x,y|t) can be expressed as: 



00 CO 



<Ku,v,0) = l l 

x=o y=o 



l pCx.ylt)© 11 ] u x V y 

t=max (x,y) 



$>(u,v,0) - s0 $(u,v,0) 



- 1 " 



1-S0 



(qu + pr + ruv 



-1 



It is to be noticed that if (x,y) £l, then the corresponding 
probabilities p and q in the above expressions do not give the required 
probabilities . 



III. PROBABILITY FUNCTIONS FOR THE POINTS ON ABSORBING BARRIERS 

A 

Probabilities P(x,y| t) and P(x,y| t) can be computed from 
formulas (1) and (2). In order to avoid early terminations of trials, 
these probabilities P(A) and P(B) should be computed and evaluated over 
sets A and B respectively, as illustrated by the following figure. 




These probabilities over sets A and B should be very close to 
each other in order to warrant a good design (i.e., trials with adequate 
data points) . 
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If P(A) is considerably larger than P(B), then an increase 
of friendly units (or decrease of the threat units) in the original 
mix is necessary to guarantee a successful trial. If P(A:) is consid- 
erably smaller than P(*B), similar policies can be adopted. 



Denote t p(x,M{t) the absorption probability with respect to 
the point (x,M)cA and denote p(N,y|t) the absorption probability 
with respect to the point (N,y)cB at the end of t-th (AT) trial. 

The expressions of these absorption probabilities are given below: 



P(x,M|t) 



1 

K(x,M,t) 



. . . im x-k M-k k t-x-M+k 
M(t-l): q p r 3 

(x-k) I (M-k)! k! (t-x-M+k) ! 



( 1 ) 



provided 0 ^ x < N and MSI. 



P(N,y|t) 



.. .x, N-k y-k k t-N-y+k 
N(t-l) ! q r s J 

(N-k)! (y-k)! k! (t-N-y+k)! 



( 2 ) 



provided 0 ^ y < M, N s 1. 



. 1N1 N-l-k M-l-k k t+l+k-N-M 

(t-1)! q p r s 

(N-l-k)! (M-l-k)! k! (t-l+k-N+2-M) ! 



provided N > 1, MSI. 



P(N,M|t) = r l 

k(N-l ,M-1 , t-1) 



IV. ABSORPTION PROBABILITIES: 

Let us denote P^(N,M) the probability that the random walk 
will terminate at the point (a,M) and by P (N,M) the probability that 

A 

it will terminate at any point on the absorbing barrier A, then we have: 

P 0 (N,M) = l P(a,M|t) 

t=max(a,M) 



provided 0 < a < N and M > 1, and 



Denoting by i(a,M) the number of updatings to be performed until the 
particle will be absorbed at point x(a,M), the distribution of x(a,M) 
should be defined by the conditional probabilities as follows: 




for t - max(a,M), max(a,M) +1, ... 



The expected number of updatings is: 




(4) 



The computation (4) is pretty combersome for the above quantity. Yet 
E[t] is very important for us to know. If it turns out that E[t], 
the expected number of updatings, is very small, we will have insufficient 
data during the trial. If E[t] is large, then the selection of (AT) 
may be too small or the trial may last too long . 

The following paragraphs will be devoted to a discussion how 
E[t] can be used in computing the expected duration of the trial. 



Denote by D^(N,M) the expected duration of the random walk, assuming that 
it will terminate at any point on the absorbing barrier A. Then we have: 



If (ct,M)eA, the quantity: 



00 



l t P(a,M|t) 



t=max(a,M) 



a M 

would be given as the coefficient of the term u v in the expansion of 
the following power series: 



1 - (qu + s) 



1 - (qu + pv + ruv + s) 



D a (N,M) = 



— l E(t(o,M)) p (N,M) . 

P (N,M) o £ 0 « 



1 
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V. APPLICATION 



Consider the rectangle illustrated in Figure 3. We wish to 
determine the mixes of forces, i.e., to determine N,M such that with 
high probability the particle will be absorbed in the absorption barrier 
in the set C* (see Figure 3) where: 

C* = {(x,y): either x = N,yiy<Mory = M, — <x<N^ , 

In other words, one can try to find various values of N,M such that the 
absorption probability has reached to some acceptable level (say 80%, 
i.e. P(C) = 80%). Once such an optimal mix has been found, experiments 
can be repeated with a change of test conditions, i.e., terrains or 
probability of kill tables, in order to obtain data which in turn will 
be used 'to detect any probability changes. As the combat between the 
friendly force and the threat force going on in the battle zone, one 
keeps a continuous record on the quantity defined as the exchange ratio 
at time t. Definition of Exchange Ratio at Time t (denote as ER(t)) is 



ER(t) = 



// of threat tank loss up to and including time t 

total # of threat tanks 

// of friendly anti-tank systems loss up to and 

include time t 

total // of friendly anti-tank systems 




(N,0) 



Figure 3 



The above definition is used for !, exchange ratio at time t n if data 
source and instruments are at high performance level, A typical graph 
of this exchange ratio ER(t) is given in Figure 4. 
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Figure 4. 



(N,M) 



( 0 , 0 ) 




. In order to control the balance of force structures and 
mixes, it is ideal to have some bounds, say a and b such that 
the probability that ER(t) is bounded in between a and b for 
all t > 0 is very high, i.e.. 



P r (a < ER(t) < b; 0 < t) = .95 

It is clear that a = tan 0^ and b = tan 0^ where 0^ and 
are shown in Figure 4. 

Using a two-dimensional random walk model [2. 4], -we have illus- 
trated above how to determine the best force structures such that the 
exchange ratio is close to one for most of the times during a trial. In 
other words, the probability that the sample path falls into the shaded 
region (see Figure 5) is maximized by an optimal choice of force structures. 



( 0 , 




Figure 5 



( 1 , 1 ) 



Due to the fact that instruments are most likely working at 
low performance levels, we provide an alternative definition for "exchange 
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ratio” using indirect measurements. This proposed definition of ”exchange 
ratio” is the ratio of expected number of kills normalized by the initial 
force units. 

The above definition of the exchange ratio is closer to the 
measure of effectiveness rather than the exchange of kill ratio but due 
to the instrument limitations this may be a good alternative measure of 
exchange ratio. 

Another way of defining the same quantity is by looking at the 
remaining kill potential or fighting value of the forces. This remaining 
fighting force would be compared with their respective fighting potentials 
relative to their respective initial fighting potential. The degradation 
of their respective potentials over time will be plotted on a graph sheet 
for each trial. 

The literature indicated that discrete state stochastic model 
formulations of combat have been difficult to solve even when the process 
is considered to be Poisson (Lanchester type) with stationary transistion 
mechanisms II]. 

The technical portion of our discussion seems new and provides 
a discrete state two-dimensional stochastic model of combat which can be 
extended into a complete stochastic model for field experimentation such 
as the Tactical Effectiveness Testing of Anti tank Missiles [3]. 
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